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ABSTRACT 


Literal expressions for the preeessional motion of the mean equator referred to 
an arbitrary epoch are constructed. Their numerical representations, based on 
numerical values recommended at the working meeting of the International Astronom- 
ical Union Commission held in Washington in September 1974, are obtained. In con- 
structing the equations of motion, the second-order secular perturbation and the secular 
perturbation due to the long-periodic terms in the motions of the moon and the sun are 
taken into account. These perturbations contribute more to the motion of the mean 
equator than does the term due to the secular perturbation of the orbital eccentricity 
of the sun. In this paper, we use the correction Af = XV 10 to the speed of Newcomb's 
lunisolar precession, the correction Ag - -0'.'03 to the speed of Newcomb's planetary 
precession, and the geodesic precession f = 1'.'915. 


RESUME 

On construit des expressions litterales du mouvement precessionnel de 
l'equateur moyen en prenant comme reference une epoque arbitraire. Leurs 
representations numeriques sont obtenues en se basant sur les valeurs nu- 
meriques recommandees aux seances de travail de la Commission de 1' Union 
Astronomique Internationale tenues a Washington en septembre 197^- On a 
tenu compte, en construisant les equations du mouvement, de la perturba- 
tion seculaire du second ordre et de la perturbation siculaire dQe aux termes 
6 longue periode dans les mouvements du soleil et de la lime. Ces perturba- 
tions contribuent plus au mouvements de l'equateur moyen que ne le fait le 
terme du a la perturbation seculaire de 1 ' excentricite de l'orbite du soleil. 
Nous employons, dans cet expose, la correction Af = 1,10" pour la vitesse de 
precession lunisolaire de Newcomb, la correction Ag =-0,03" pour la vitesse de 
la precession planitaire de Newcomb et la precession geodesique fg = 1,915" 


v 



KOHCTIEKT 


CocTaBJieHH SyKBeHHbie BwpasceHUH npeneccHOHHoro flBHaceHMH 
CpeflHerO 3 KBSTOpa OTHOCHTe JIBHO np0K3B0JIBH0 3anaHHOM anoxH. 
IIojiyueHH mx tiMCJieHH&ie BbipaxcsHHH , ocHOBaHHbie Ha qwcjieHHbix 3Ha- 
ueHHHX, KOTopbie (5mjim peKOM eHflOBaHb! Ha patfoueM coBemaHMH kommccmm 
MescflyHapoflHoro acTpOHOMHnecKoro cotoaa, coctohbih6mch b ceHTH<5pe 
1974 r. b BamMHTTOHe. ripw cocTaBJieHHM ypaBHeHHH flBwxeHHH 6 hjim 
ytiTeHbi BeKOBoe B03MymeHHe BToporo nopnziKa m BeKOBoe BOSMymeHHe , 
BbI3BaHH0e HOJironepHO&HtjeCKMMM UJieHaMH B ypaBHGHMflX flBMXeHHH 
nyHbi m cojiHua. 3tm bo 3MymeHMH <5ojiBiue bjimhjot Ha ypaBHenne flBrae- 
hhh cpeHHero 3 KBaTopa , yeM yjien, ynHTbtBaiomHM BeKOBoe B03MymeHHe 
op^HTaJiBHoro 3 KcneHTpHCHTexa cojiHna. B HacTonmen paCoTe mh nojiB- 
3yeMCH nonpaBKOM Ha ckopoctb uyHHo-cojiHeyHOH npenecCHM HBK>KOM<5a 
Af = 1,10", nonpaBKOM na ckopoctb njiaHeTapHOM npeneccMH HBWKOMtfa 
Ag = -0,03" m reone3HnecKOM npeneccHefi f g = 1,915". 


vi 



FORMULAS FOR PRECESSION 
Hiroshi Kinoshita 

I. HAMILTONIAN FOR PRECESSION 

In this paper, we basically adopt the conventions and notations used by Andoyer 
(1911). The Hamiltonian K for the precessional motion of the mean equator is as 
follows (Kinoshita, 1975): 

K 
F 

s 


E 

where H = G cos I and G= C«. The angle J between the axes of the figure and the 

6 

angular-momentum axis of the earth is of the order of 10 , Therefore, we can neglect 

3 2 

the term — sin J in this paper* Here, the ecliptic at initial epoch t = 1850 is adopted 

a 

as the reference plane of the inertial system, and the mean equinox at 1850 is the 
origin for longitude. The longitude h of the mean equinox of date is measured from the 
fixed mean equinox of 1850 along the fixed ecliptic of 1850 and then along the ecliptic 

This work was supported in part by Grant NGR 09-015-002 from the National 
Aeronautics and Space Administration. 
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of date, and the inclination I of the mean equator is referred to the ecliptic of date 
(see Figure la); h and I, respectively, correspond to the general precession in longi- 
tude and the obliquity, but their signs are opposite from conventional usage, since we 
are using a right-hand system here. The principal moments of inertia of the earth 
are denoted by A, B, and C, and the angular velocity of the rotation of the earth is w. 
The masses of the sun, the earth, and the moon are denoted by Mq, M^, and M^; the 
sidereal mean motions of the sun and the moon by n^ and ; and the mean motion of 
the node of the moon by n . The first term F in the Hamiltonian K is the secular term 

due to the direct lunisolar perturbation. The factors M ft , Mi, and M 0 in the equation 

13 9 u * 

for F are the secular terms of -r (a/r) (1-3 sin p), (a/r)° sin p cos p sin X., and 

S o & 

(a/r) 3 cos J p cos 2\, respectively, where X and p are the longitude and latitude of the 
moon referred to the ecliptic of date and the mean equinox of date. The factors and 
Sg are the same quantities as for the sun, and AS^ is due to the secular perturbation of 
the orbital eccentricity of the sun. The term with Mg comes from the second-order 
secular perturbation; The quantities M^, M^, Mg, Mg, Sq, Sg, and AS^ depend only 
on the orbital elements of the moon and the sun, and their numerical values (Kinoshita, 
1975) are obtained from Brown's theory of the moon as improved by Eckert, Walker, 
and Eckert (1966) and Newcomb's theory of the sun. The terms having M^, Mg, and 
Mg as factors are not included in Newcomb's preeessional theory; Mj, Mg, and Sg come 
from the long-periodic terms in the motions of the moon and the sun. The second term, 
E, in the Hamiltonian K represents the effect due to the motion of the ecliptic, which 
is caused by planetary perturbations. The form of F is independent of whether or not 
the reference system is inertial. Newcomb (1906) gave the motion of the ecliptic of 
date in the following form: 

sin k sin <() = pt + p't 2 + p"t 3 - 57341t + 071935t - 0’.'000l9t , 

( 2 ) 

sin k cos <f> = qt + q't 2 + q"t 3 = -46’.’838t + 0'.’563t + 0700035 , 

in which t = tropical centuries from 1850, k is the inclination of the ecliptic of date on 
the ecliptic of 1850, and <p is the geocentric longitude of the heliocentric ascending node 
of the ecliptic of date of 1850. If the ecliptic does not move, the additive term E dis- 
appears, and both the obliquity and the rate of the general precession in longitude 
are constant because the canonical variable h conjugate to H is cyclic. 
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2. EQUATIONS OF MOTION AND THEER SOLUTIONS 


Equations of motion are given from the Hamiltonian (1) as 


dh _ 9K dH _ 9K 
dt 3H ’ dt ~ 9h » 

or (3) 

dh 1 9K dl 1 9K 

dt Cw sin I 91 ’ dt Cu> sin I 9h 

Expressing the mean obliquity as £ (£ = -I) and the general precession in longitude as 
(^= -h), which is reckoned westward, equations (3) become 


(if = shi^ £ " cot £ sin + 4>) ^ + sin k cos + <j>) - (1 " cos k) 


d<t> 

dt * 


^ = cos + 4>) "jj~ - sin k sin + 4*) ^ 

= cos \f> ~ (sin k cos <f>) - sin * (sin k sin 4>) + 2 sin 2 cos + 'l 5 ) ^ > 


R(£) = (Pq + P^t) sin £ cos £ + cos 2£ + Vq( 6 cos £ - 1) sin £ , 


p o= k «( M o-I M 2 ) +k o( s o-| s 2 ) 


Qo-\ M i > 


PRECEDING PAGE BLANK NOT FILMED! 
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Although equations (4) and (5) are nonlinear, their solutions can be obtained by a 

successive approximation, assuming that 4* and £ can be expanded into power series of 

time with the initial conditions = 0 and £ = £. at t = 1850. After we solve for and 

0 

£, the other quantities — the obliquity referred to the fixed ecliptic £ ^, the lunisolar 
precession in longitude ip, and the planetary precession x — are derived with the aid 
of the trigonometric relation in the triangle y^Ny (see Figure lb). To simplify the 
expressions and for convenience in systematic numerical calculations, we solve the 
following equations together with equation (5): 


d£ 

= R(£) sin x , (6) 


sin 


£ 1 Ht = cos x ’ 


( 7 ) 


sin £^ sin x - sin k sin (\j> - cf>) , 


( 8 ) 


tan | (* - $) 


- tan x COS (£ + t^/2 
cos (£ - S^/2 


(9) 


where equations (6) and (7) are easily derived from equations (4) and (5) with the aid 
of the trigonometric relation in the triangle y^N-y. Equations (6) and (7) represent the 
motion of the pole of the mean equator in the ecliptic system and are equivalent to 
those obtained by Andoyer (1911) if R(£) is replaced by P cos £, or, in another words, 
ifQ 0 .Vo = °. 
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We can express £, £^, ip, x> and R(£), la the form 

£ = £q + at +• a't 2 + a"t 3 , 

£ 1 = £ 0 + bt + b't 2 +b"t 3 , 

ip = ft 4- f't 2 + f"t 3 , (10) 

X = gt + g't 2 + g"t 3 , 

* = ht + h't 2 + h"t 3 , 

R(£) = R 0 (£ 0 ) + + R 2 (€ 0 )t 2 . 

Substituting these into equations (5) through (9), and equating the coefficients of like 
powers of t, we get, to the third power of t, 

a - q j 

a' = q'- |ph , 

a" = q" - ~ q(h 2 - p 2 - q 2 ) - | (2p'h + ph') , 


b = 0 , 


b ' = l R o g = i pf » 

b" = | (R^ + R Q g') = I (g'f + 2gf') sin E q 


( 11 ) 


g = p cosec £ Q , 


g' = (p' + qh) cosec £ Q , 

g" = (p" -t- qh' 4- q'h - ~ ph 2 ) cosec £ Q + £ g 3 - gb' cot £ Q , 


h = f - g cos £ Q , 


[eq. (11) cont. on next page] 
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h' = f' - g' cos £ Q + 2 pq 5 

h" = f" - g" cos £ q + g(a' + b') sin £ 0 + ^ ag' sin £ Q - p 3 cot £ Q , 


f = R q cosec £ Q 

2 

= P 0 cos £ 0 + Qq cos 2£ 0 cosec £ Q + V Q (6 cos £ Q - 1) , 

f' = cosec £ Q 

= ^ P 0 a cos 2£ Q cosec £ Q + f cos £q _ 2 % a cos ^0 ( n ) 

+ V 0 a{9 cos 2 £ q - 4) cot £ 0 , 

f" = | (R 2 cosec £ q - | fg 2 - fb' cot £ Q ) 

= “ |p 0 {a' cos 2 £ q - a 2 sin 2£ Q ) + P^ cos 2£ Q - 2Q Q (a' sin 2£ 0 +■ a 2 cos 2 £q) 

+ V^a' ( 9 cos 2 £ q “ 4) cos £ Q - ^-a 2 sin £q (14 + 27 cos 2£ Q )J | cosec £ Q 
-ifg 2 4fb' cot £q , 

which can be written in many different ways. We have chosen these in such a way that 
the quantities Pq, P^, Q q , and V Q appear only in the lunisolar precession in longitude. 

As a result, except for the lunisolar precession, they coincide with Andoyer’s (1911). 

It is worth noting that the obliquity referred to the fixed ecliptic does not have a 

2 3 

secular term in the first power of time, although it does have the terms t , t , « • • . 
However, Woolard (1953) derived a secular term -0i’0256t, which has an effect amount- 
ing to 10% of the generally accepted value of Ort’s constant B. This discrepancy occurs 
because Woolard neglected some terms in constructing the equations of motion 
(Kinoshita, 1975). 
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According to Newcomb's measure (1906), the general precession in longitude is 
the longitude of the mean equinox of date, referred to the fixed mean equinox and 
ecliptic of epoch: <J< N = i p - y^L, where L is the orthogonal projection of the mean 
equinox of date onto the fixed ecliptic (Figure lb). This measure is favored by Plummer 
(1916) and Woolard and Clemence (1966). From the triangle y^Ly, we have 

tan'\~H = tan x cos £^ (12) 

and then, to the third order of t, 

% " V + h k ‘ 2 + h &* 3 ’ < 13) 


where 


h N = f - g cos £ q , 

hjj “ f ' - g ' cos £ 0 ’ ’ (14) 

hj = f" - g" cos £ q + gb' sin £ q - ± p 3 cot £ Q . 

The expressions for the general precession in longitude according to Newcomb's meas- 
ure differ in the second and third orders of t from those [see eqs. (11)] according to 
Andoyer's definition, and the numerical differences between them are only of the order 
of OV0005t 2 and O’.’OOOOlt 3 . 

Formulas for <f> and k can be obtained from equations (2): 

4 > = 4> 0 + 

k = wt + 


vt + v't 2 


w't 2 + w"t 3 , 


(15) 
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where 


*0 = tan q > 

p ' cos <t> 0 - q' sin <t> Q 


p" cos - q" sin <f> 0 - w'v 

V' w ' 

(16) 

w = Vp 2 + Q 2 > 


w' = p' sin <j> 0 + q' cos <j> Q , 

1 2 13 

w" = p" sin (j)^ + q" cos 4> 0 + - wv 

Next we calculate the precessional quantities in the equatorial system. From 
equations (4) and (5) and with the aid of the trigonometric relation in the triangle 
(Figure lb), the motion of the pole of the mean equator in this system is 
expressed as 

sin j ^ = R(£) cos p , (17) 


= R(£) sin p . 


(18) 


The spherical triangle (Figure lb) gives the following relations: 

u + p+v -tan (*/2) sin[(£ 1 + £ 0 )/2] 
tan 2 “ sin \WP £ q )/2] ’ 


(19) 


tan 


tan (0/2) cos [(E 1 + £ Q )/2] 
cos [(£j - £ 0 )/2] 


( 20 ) 


As t approaches the epoch in equations (19) and (20), we have tan (p + p)/2 — <*> and 
tan (p - p)/2 — 0, because x — gt, ip — ft, (£j + £ Q )/2 — £ Q , and (£ 1 - £ Q )/2 — b't 2 /2. 
Therefore, we can express p + p, p - p, and j as follows: 
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( 21 ) 


p + p = 180 + ut + u't 2 + u"t 3 , 
p - p = rt + r't 2 + r"t 2 , 
j = st + s't 2 + s"t 3 . 

Substituting these into equations (17), (18), and (20), and equating the coefficients of 
like powers of t, we get, to the third power, 


s = R q = f sin £ Q , 
s' =|R 1 = f'sin £ q , 

s" = | r 2 " 24 R o r2 = sin £ 0 + 24 fshl £ 0 ^ + 6fg C0S £ 0 " ^ COs2 £ 0 J 5 


u = 0 


sr' - rs' 
3s 


sr" - 2u' s' - rs" _ rs r 
2s 24 48 


( 22 ) 


r = f cos £ Q - g , 
r' = f' cos £ Q - g' , 

r" = f" cos £q - s" + ^ ^ sin2 ^0 cos £ 0 " 3 S) • 

The quantities £, and z introduced by Newcomb (1906) are related to p, 
p, and xby^-p + x~ 90% = 90 ° ~ P ’ and z = M “ 90 °' 

In order to obtain formulas for epoch t' referred to an arbitrary epoch t, we must 
first derive the position of the ecliptic of epoch t' referred to the ecliptic and the mean 
equinox of t in terms of the coefficients for the fundamental epoch t = 1850. Denoting 
the quantities for epoch t' referred to an arbitrary epoch t by (k), (<b), (£), • • ■ and 
their coefficients by (p), (p'), (p"), (q), (q'), (q"), (a), (a'), (a"), we have 
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sin (k) sin ((j)) = (p) 9 + (p') 0' 2 + (p") 0' 3 , 

sin (k) cos (<J>) = (q) 0 + (q') 0 2 + (q") 0 3 , 

(£) = (E q ) + (a) 0 + (a') 6 2 + {a") 0 3 , (23) 

« 

* ? 

where 0 = t' - t. It is clear that (P Q ) = P Q + P^ and (£ Q ) = £. From the triangle 
N(N)N' (Figure 2), we get 

(p) =p + Pl t + p 2 t 2 , (q) =q + q A t + q 2 t 2 , 

<p'> = p' + p' A t , (q'l = q' + q' x t , (24) 

(P") = P" , <q")=q" , 

in which 

Pj = 2p' + qh , q x = 2q' - ph , 

p 2 = 3p" + 2q'h + qh' - | p(h 2 - p 2 - q 2 ) , q 2 = 3q" - 2p'h - ph' - 1 q(h 2 - p 2 - q 2 ) , 

pi = 3p" + q'h + |p(p 2 + q 2 ) , q' : = 3q" - p'h + |q(p 2 + q 2 ) , 

(25) 

which were derived by Andoyer (1911). 

The relationships among the coefficients for epoch t' referred to an epoch t are 
the same as those for t referred to the fundamental epoch t = 1850: 
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411 — 09 J 

















(a) = a + a^t + a^t 2 , 

(a') = a' + a'jt , 

11 


g 

tl 

O 

V* 

<b') = b' + b^t , 

•v' 

II 

cr 

(26) 

(g) = g +■ g^ + g 2 t 2 , 

(g') = g' + gjt , 

(g") = g" , 



• } 

where 



a l = q l " i (ph l + p i h) » 

b l = 1 (Pl f + pf l ) » 

2 

g^^ = - ap cos £ Q cosec £ Q + p 1 cosec £q , ^ 

g 2 = 2 P [ (2 “ ^ £ 0 ) ^ " a ' sin2 € o] cosec3 £ 0 " p I C0S £ 0 cose ° 2 £ o 

+ P 2 cosec £ Q , 

g' x = -a(p' + qh) cos £ 0 cosec 2 £ Q + (p^ + qh^ + q^h) cosec £ Q , 

- g y cos £ q + ag sin £ Q , 

h 2 = f 2 “ & 2 COS £ 0 + g(a ' Sin £ 0 + 2 ^ C0S £ 0 ) + S 1 a Sin £ 0 ' 

h i = f i “ g i cos £ o + g ' a sin £ 0 + I ^* q l + p l^ » 

[eq. (27) cant, on next page] 
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*1 = q ^1 C ° s2 *0 " q i Sin ^0 COS V ’ 

4>2 = ^ (p 2 cos2 *0 “ 4>i<i tan <t> 0 " q i^i " q 2 sin cos V » 

V 1 = w [ P 1 C ° S +0 " p,< ^l Sln +0 ” q ^l C0S ^0 “ q i Sin *0 


IT (p ' cos ' 


0 - q ' sin V] ’ 


W l = w* >p l + qq l ) ’ 

1 (A + < h \ 1 2 

w 2 = w \—T- + PP 2 + q V - ^3 (P p l + qq i> ’ 

w' x = p'$ 1 cos <|> 0 + p' A sin <t> 0 - q'c^ sin <j> Q + q' cos , 


r = f JL cos £ Q - af sin £ Q - ^ , 


r 2 = -f (a' sin £ Q + ^ a 2 cos £ Q ) - a£ 1 sin £ fl - g 2 + f 2 cos £ Q , 
r j = f' x cos £ Q - f'a sin £ Q - g^ , 


Sj = f^ sin £ q + fa cos £ Q , 


1 2 

s 2 = f 2 sin £ 0 + f^a cos £ Q + f(a' cos £ Q ~ j a sin £ Q ) , 


s'j = f^ sin £ q + f'a cos £ Q , 

u i = i (sr i * s i r ' - s ' r i ' s i r ' 3 u ' s i> > 

1 2 

f^ = Pj cos £ 0 ~ P Q a sin £ Q + - Q Q a (cos 3£ Q - 3 cos £q) cosec £ Q 


6V q a sin 2£ 0 , 


[eq. (27) cont. on next page] 
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f 2 = - P Q (a' sin £ 0 + 1 a 2 cos £ Q ) - sin £ Q - 6V Q (a' sin 2£ Q + a 2 cos 2£ Q ) 

I" 2 3 „/ 

+ Q j(9-2 cos 2£ 0 + cos 4£ Q ) cosec £ Q -g- * cos 3£ 0 “ 3 cos ^ 

2 c 1 

cosec £ 0 , 

J (27) 

f' x = - I P 1 a sin £ 0 - 1 P Q aq (4 cos £ Q + cos 2£ Q cos £ Q cosec 2 £ Q ) 

+ 1 (P x q + P^) cos 2 £ q cosec £ Q - 2 Q q ^ cos £ Q - q 2 sin £ Q ) 

+ 1 V 0 a(9 cos 4 £ q - 18 cos 2£ Q - 1) cosec 2 £ Q . 
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3. NUMERICAL REPRESENTATION OF PRECESSIONAL MOTION 


In order to obtain the numerical representation of precessional motion, it is 
necessary to determine values of P^, P^, Q^, and V^, which depend on the value 
of the speed of the lunisolar precession. We adopt the following: 


f 0 = 5037'.'08 + IV 10 + 1V9I5 = 5 040 '.’09 5/tropical century, 

M_ - = 496303. 4 X 10~ 6 , 

U Z J 

M = -20. 7 X 10~ 6 , 

M = 2992X 10 _6 , 

d 

(S 0-1 S 2 ) 1900 = 500210 - 0X J °‘ 6 > 

AS Q = -1.056 X 10“ 6 , 

M./M =0.02130 , 

Mq/(M 0 + M^) = 328900. 2 , 

n^/w = 6 3 238 90 8'.’ 5/tropical century, 

n^/w = 35386'.'791/tropical century, 

n^/wn Q = -9. 082445 , 

F = 0. 9990931420 , 

8 1900 = 23 ° 27 ’ 8 " 28 » 

Ap = -0'.'012 , Aq = Oi'037 . 


(28) 
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The rate of lunisolar precession at 1900, f 0 , is obtained by applying the proper- 
motion correction Af = IV 10 (Fricke, 1971) and the geodesic precession f = IV 9 15 
(de Sitter, 1938) to Newcomb’s value 5037V08. Values of M Q , Mj, M 2 , Mg, S 0’ S 2’ 
and AS Q were calculated by Kinoshita (1975) by using Brown’s theory of the moon 

as improved by Eckert et al . (1966) and Newcomb's theory of the sun. Values for 
/M 0 and Mq/(M^ + Mg) were recommended at the working meeting of the Inter- 
national Astronomical Union held in Washington in September 1974. The meeting 
also recommended a new system of planetary masses. Since the motion of the ecliptic 
depends on the planetary masses, a change in mass causes a change in p, q, p', q', 
p", and q" in equations (2). We adopt here Ap= -0V012 and Aq = 0V037 as corrections 
to Newcomb's p = 5V341 and q = -46V838, which were calculated by Laubscher (1972) 
from the planetary masses proposed by Klepczynski, Seidelman, and Duncombe (1971). 
Klepczynski et al . 's values are slightly different from those recommended, but the 
changes are too small to affect Ap and Aq. The correction to planetary precession 
due to Ap is -0V03 per tropical century. The value of the obliquity includes the 

effect of Aq. Using values from equations (28) and solving the following equation 
[see eqs. (11)] with respect to the dynamical ellipticity of the earth [2(-A - B)/2C], 


’woo- 5040 - 095 


we get 


= 3X 


2C - A - B 
2C 


M„ 


. _L JL 

M { + M 0 F 9 “ 


( M o - 5 m2 ) 


cos £ 




1VL 




X 2C - A - B ^ _2 ^ „ 

, ~ 1 ” — ~ (6 COS C “ 1) 

1 sin £ 4 IVL + 2C 


M, 


n. 




2C - A - B 
2C 


= 0.003273952 , 
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and 


kj. = 7567’.’55998/tropical century, 
k Q = 3475V32347/tropical century, 

P g = 5494’.' 1988/tropical century, 

(29) 

= -O'.'00367/tropical century 2 , 

Q g = -0'.' 1566/tropical century, 

Vg = 01'008 18/tropical century. 

We see that Q g and V Q are larger than P^t within the period in which Newcomb's theory 
for the ecliptic is valid. Strictly speaking, as seen from equations (10) and (11), the 
obliquity at 1900 depends somewhat on P g , Pp Q g , and V Q through the second and 
third orders of t, but the effect on £ due to them amounts to only 01'002 and, 
therefore, does not have a numerical effect on the determination of the dynamical 
ellipticity and of P Q , P , Q q , and V Q . 

Applying equations (29) to equations (11), (16), (22), (25), and (27), we derive 

numerical representations for precessional motion; these are given in Table 1. 

Table 2 shows the precessional formulas, in which t' and t are reckoned from 2000 
d h 

January l a 12 ET in units of Julian centuries. The geodesic precession which is 
a direct motion of the zero point of longitude along the ecliptic of date, is subtracted 
from f, h, r, and s in equations (11) and (22) to obtain f g , h g , r g , and s g in Tables 1 
and 2, as follows: f Q = f - f ^ h Q = h - f^ r Q = r - f g cos £ g , and s g = s - f^ sin £ g . 
For comparison, Tables 3 and 4 show precessional formulas based on Newcomb's 
values of P Q = 5490V66, P 1 - -0'.’00364, and Q q = V Q = 0. 
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Table 1. Precessional formulas: 6 = t' - t, where t' and t are reckoned from 1850 in units of tropical centuries. 


sin (k) sin (<j>) = (5'.’3290 - 0’.'75374t + 0».'000335t 2 ) 0 + (0'.'19350 + 0’. T 000802t) 0 2 - O'.'OOO19O0 3 

sin (k) cos <c|>) = (-46 '.'80 10 - 0'.'01729t + 0'.’005489t 2 ) 0 + (0'.'05630 - 0'.'003668t) 0 2 + O'.'OOO35O0 3 

(<}>) = 173°30 T 14'.'36 + 3288"00125t + 0'.'576239t 2 + (-869'.’78374 - 0’’554964t) 0 + O'.'O245760 2 

(k) = (47'.' 1034 - 0'.'06809t + 01'000569t 2 ) 9 + (-0’.'03405 + 0’.'000569t) 0 2 + 0’. T 0000509 3 

(£) = 23°27'31'.'68 - 46'.'80100t - 0’.’008645t 2 + 0'.'00183t 3 + (-46'.'8010 - 0l’01729t + 0'.'005489t 2 ) 0 
+ <-0i'00865 + O’. , OO549Ot)0 2 + O'.'OO183O0 3 

(£^) = 23°27 T 31V68 - 46'.'80100t - 0'.'008645t 2 + 0l'00183t 3 + (0'. T 06510 - 0’.'009202t) 0 2 - OVOO77330 3 
ty) = (5037’.'9337 + 0'.'49264t - 0'.'000039t 2 ) 0 + (-11'05099 - 0”001513t) 0 2 - O1'OO153O0 3 
© (x) = (13'. , 3864 - l'.'88640t - 01’000144t 2 ) 0 + (-21'37947 - 01'001554t) 0 2 - Oi'OO16610 3 

(*) = (5025'.'6537 + 2 '.' 22191 t + 0'.'000264t 2 ) 0 + {r.'13120 + 0 T .'000212t) 0 2 + O’. , OOO1O2 0 3 
(p) + (p) = 180° + (0”79223 + 0'.’000648t) 0 2 + O'.'OOO3280 3 

(p) _ ( p ) = (4608'.' 1448 + 2'.'79355t + 0'.'000118t 2 ) 0 + (1'.'41535 + 0'.'000071t) 0 2 + O'.'O363930 3 
(j) = (2005'.'5487 - 0’.'85290t - 0’.'000363t 2 ) 0 + (-0'.'41839 - 01'000384t) 0 2 - O'.'O41877 0 3 



Table 2. Precessional formulas: 0 = t' - t, where t' and t are reckoned from 2000 in units of Julian centuries. 


sin (k) sin (<fy = {4’.'1992 - 0'.’75277t + 0'.’000335t 2 ) 0 + (Of 19471 + 01'000802t) 0 2 - OL'000190 0 3 

sin (k) cos (<j>) = (-4018156 - 0'.'00082t + 0f005489t 2 ) 0 + (0'. f 05080 - 0'.’003668t) 0 2 + 0”0003500 3 

(<(>) = 174°52 , 27'.'66 + 3289L'80023t + 0l’576264t 2 + (-8701'63478 - 01'554988t) 0 + O’.'O245780 2 

(k) = (47V0036 - 0'.'06639t + 0'.'000569t 2 ) 9 + (-OL'03320 + 01'000569t) 9 2 + O”OOOO5O0 3 

(£) = 23°26 , 21’. , 47 - 46'.’81559t - 0:'000412t 2 + 0'’00183t 3 + (-461'8156 - 0L’00082t + 01’005489t 2 ) 0 
+ (-OV00041 + 0”005490t) 0 2 + O1'OO183O03 

(£^) = 23°26 , 211’47 - 46'.'81559t - 0l’000412t 2 + 0'.'00l83t 3 + (0'/05130 - 01'009203t) 0 2 - Ol f OO77340 3 

tyj) = (5038!' 7802 + 0 f .'49254t - 0'.'000039t 2 ) 0 + <-ll T 05331 - 0l’001513t) 0 2 - O’’OO153O0 3 

(X) = (10”5567 - ll'88692t - OV000144t 2 ) 0 + (-2'.’38191 - 0'’001554t) 0 2 - O'.’OO16610 3 

(*) = (5029”0946 + 2i'22280t + 0V000264t 2 ) 9 + (1'.' 13 157 + 0’/000212t) 0 2 + O'.'OOO1O2 0 3 

(|a) + (p) = 180° + (0’.'79323 + 01’000648t) 0 2 + Oi'OOO328 0 3 

(p) _ ( p ) = (4612V4339 + 21’79402t + 0’.'000118t 2 ) 0 + (l'. r 41551 + 0'.’000071t) 0 2 + O’.'O363950 3 
(j) = (2004"3113 - 01'85403t - 0'.'000363t 2 ) 0 + (-0'. r 41898 - 01'000384t) 0 2 - O'.'O4188O0 3 



Table 3. Precessional formulas based on Newcomb's values; 9 — t r - t, where t' and t are reckoned from 1850 
in units of tropical centuries. 

Sin (k) sin {$) = (5'.’3410 - 0f75393t + 0f000336t 2 ) 0 + (Of 19350 + 0f000802t)9 2 - 0f0001909 3 

sin (k) cos ( 4 .) = (-46f8370 - 0f01750t + 0l’005489t 2 ) 0 + (Of 05630 - 0f003665t)9 2 + 0f0003509 3 

(<()) = 173°29'40f01 + 3286f28345t + 0f560350t 2 + (~869f 12350 - 0f550988t) 9 + Of 025375 0 2 

(k) = (47 f 1405 - 0f06803t + 0f000567t 2 )9 + (-Of 03401 4- 0f000567t)9 2 + Of 000050 0 3 

(£) = 23 e 27'31f68 - 46f83700t - 0f008752t 2 + 0f00183t 3 + (-46f8370 - 0f01750t + 0f005489t 2 )9 
+ (-Of 00875 + 0f005489t) 9 2 + 0f00l830 9 3 

(£ x ) = 23°27’31f68 - 46f83700t - 0f008752t 2 + 0f00183t 3 + (0f06521 - 0f009199t)9 2 - 0f0077309 3 

(ip) = (5036f 8381 + 0'.'49299t - Of 000037t 2 ) 9 + (-lf07129 - OfOO1478t)0 2 - OfOO15330 3 

(X) = (13 '.'4166 - lf88685t - OfOOO142t 2 )0 + (-2f37995 - Of 00157 It) 9 2 - Of 001657 9 3 

( 4 -) = (5024f5305 + 2f22267t + Of 000264t 2 ) 9 + (If 11134 + Of 000263t) 9 2 + 0f0000959 3 

<p) + (p) - 180° + (Of 79236 + 0f000656t) 0 2 + Of 000328 0 3 

(|j) _ (p) = (4607f 1096 + 2f79440t + 0f000118t 2 )9 + (lf39720 + 0f000118t)9 2 + 0f0363209 3 
0) = (2005f 1125 - Of 85294t - 0f000365t 2 ) 0 + (-0f42647 - 0f000365t) 0 2 - OfO418O2 0 3 



Table 4. Precessional formulas based on Newcomb’s values: 0 = t' - t, where t f and t are reckoned from 
2000 In units of Julian centuries. 



to 

CO 


sin (k) sin (<f>) = (41*2109 - 01’75296t + 01'000336t 2 ) 9 + (01*19471 + 0*’000802t) 0 2 - 01’0001909 3 

sin (k) cos (<f>) = (-461*8519 - Ol’OOKMt + 01’005490t 2 ) 9 + (01*05080 - 01'003665t) 0 2 + 01*000350 0 3 

(<|>) = 174°51’50l*70 + 32881*03472t + 01’560374t 2 + (-8691*96857 - 01*5510111)9 + O1*O253760 2 

(k) = (471*0408 - 01*066331 + 0 r .'000567t 2 ) 0 + (-01*03316 + 01*000567t) 0 2 + 01'000050 9 3 

(£) = 23°26'211’41 - 461'85191t - 0"000519t 2 + 0’. T 00183t 3 + (-46'.'8519 - 01*00 104t + 01*005490t 2 ) 0 
+ (—01*00052 + 0’.'005490t) 0 2 + 01*001830 0 3 

= 23°26*2U*41 - 46'.'85191t - 01*000519t 2 + 0’.*00183t 3 + (01*05142 - 0'.'009199t) 9 2 - O'.'OO77310 3 
(ip) = (50371*6851 + 0!*49290t - 0l'000037t 2 ) 0 + (-11*07355 - 01*001478t) 0 2 - O1'OO15330 3 
(X) = (101*5862 - 11*8873 6t - 01'000142t 2 ) 0 + (-2**38240 - 01'001571t) 0 2 - 0**001657 9 3 
(4<) = (5027'.’9724 + 2'.'22356t + 01*000264t 2 ) 0 + (11*11178 + 01*0002631) 0 2 + 0’.*000095Q 3 
([i) + (p) = 180° + (0**79338 + 0*.*000656t) 0 2 + 01*000328 0 3 

(|J) - (P) = (46111*3999 + 2*.’79488t + 01*000 118t 2 ) 0 + (1**39744 + 01'000118t) 0 2 + Ol'O363220 3 
(j) = (20031*8751 - 0**85407t - 01*000365t 2 ) 0 + (-01*42704 - 01’000365t) 0 2 - O1'O418O50 3 
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